Phase transition in a triplet process 
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We argue that the reaction-diffusion process 3A — > 4A, 3A — + 2A exhibits a different type of 
continuous phase transition from an active into an absorbing phase. Because of the upper critical 
dimension d c > 4/3 we expect the phase transition in 1+1 dimensions to be characterized by 
nontrivial fluctuation effects. 
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The classification of continuous phase transitions far 
from thermal equilibrium is one of the most challenging 
tasks of modern statistical physics 0. Within this field 
many studies are concerned with phase transitions from 
a fluctuating active phase into one or several nonfluctuat- 
ing absorbing states, which are believed to be associated 
with a finite number of universality classes || . For such a 
phase transition to occur it is necessary that (a) at least 
one absorbing state is dynamically accessible, (b) there 
are two competing processes for particle creation and re- 
moval, and (c) there is a mechanism which prevents the 
particle density from diverging. 

In many cases, the critical behavior close to the tran- 
sition is characterized by simple power law scaling. In 
sufficiently high dimensions the critical exponents are 
given by their mean-field values, whereas below a cer- 
tain upper critical dimension d c fluctuation corrections 
have to be taken into account, leading to nontrivial ex- 
ponents and scaling functions. For this reason the study 
of fluctuation effects in low dimensional, especially 1+1- 
dimensional systems is particularly interesting. 

The most important universality class of absorbing 
phase transitions is directed percolation (DP), which 
occurs in all processes following the reaction-diffusion 
scheme A «-> 2A, A — » 0. The critical exponents, es- 
pecially in one spatial dimension, are known to a very 
high precision |j. The critical behavior of DP can be 
described in terms of a renormalizible field theory which 
was originally introduced in the context of high energy 
physics Q]. DP is relevant for experimental applications 
such as catalytic reactions || , flowing sand 0, and spa- 
tiotemporal intermittency of magnetic fluids ^ . 

The other established class is the parity-conserving 
(PC) universality class. This type of critical behavior 
is observed in a large variety of models which can be 
divided into two groups. The first group includes all 
parity-conserving particle processes || ||, [To], |ll| such as 
branching-annihilating random walks with two offspring 
A — > 3A, 2A 0. The second group of models comprises 
spreading processes with two symmetric absorbing states, 
including kinetic Ising models [fL2[ , interacting monomer- 
dimer models as well as generalized versions of the 
Domany-Kinzel model and the contact process p3. In 



higher dimensions the second group of models describes 
branching-annihilating interfaces and can be associated 
with the voter universality class |lq| . Only in 1+1 di- 
mensions the two classes of models exhibit the same type 
of critical behavior. 

Recently the pair contact process with diffusion 
(PCPD), also called annihilation-fission process, at- 
tracted considerable attention. The PCPD is a binary 
spreading process following the reaction-diffusion scheme 



nA^(n+l)A, nA — > mi, 



(1) 



with n — 2, m < 1. It exhibits a continuous phase tran- 
sition and thus could serve as a candidate for another 
independent universality class. The PCPD was already 
suggested in 1982 by Grassberger |[(|, but it took almost 
20 years until Howard and Tauber presented a first sys- 
tematic study of a bosonic variant of the process [tf7|. 
Using field-theoretic methods they were able to prove 
the existence of a phase transition, although the corre- 
sponding field theory turned out to be unrenormalizable. 
More recently, several authors studied various fcrmionic 
variants of the PCPD || [L9[ || §| || ||, |§. 
Meanwhile there is a general consensus that the critical 
behavior of the PCPD is different from all other pre- 
viously known universality classes. However, it turned 
out to be extremely difficult to estimate the critical ex- 
ponents in a reliable way, mainly because of unusually 
strong deviations from ordinary power-law scaling. 

In the present study we investigate the question 
whether further different types of critical behavior will 
emerge for n > 2. In particular we will focus on the case 
n = 3, called triplet process. As will be shown below, 
we argue that this process exhibits yet another different 
type of critical behavior. 

a. Mean field approximation. In order to determine 
critical dimension d c and the mean field critical expo- 
nents of the reaction-diffusion process ([!]) for general 
m < n, let us consider a simple mean field theory. We 
expect this process to be described by the Langevin equa- 
tion 

cV(x, t) = ap"(x, t) - p n+1 (x, t) + DV 2 p(x, t) + C(x, t) , 

(2) 
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FIG. 1: One-dimensional triplet process at criticality starting 
with a fully occupied initial state. 



which for n = 1 reduces to the well-known Langevin 
The first term accounts for both 



equation for DP 26 



particle creation and removal so that the parameter a 
plays the role of the reduced spreading probability p — p c . 
The second term is the most relevant contribution pre- 
venting the particle density from going to infinity, while 
the third term describes nearest-neighbor diffusion. 

The noise £(x, f ) takes the stochastic nature of particle 
creation and removal into account. Its amplitude has to 
depend on the local density p(x, t) since in the absorbing 
state p = there are no density fluctuations. Thus it is 
near at hand to expect noise correlations of the form 

(c(x, t)c(x', o> = iy*(x, t)s d (* - x')*(« - o (3) 

with an unknown exponent ji. For DP (n = 1), where 
the squared noise amplitude is proportional to the den- 
sity of particles, this exponent is given by fi — 1. For 
n > 1, however, the situation is more involved. Without 
the branching process, i.e., deep in the inactive phase, 
the squared noise amplitude is proportional to 9*p(x, t), 
hence p = n. At the transition, however, the branching 
process may lead to positive correlations among the par- 
ticles, increasing the intensity of the noise and thereby 
reducing the value of p. At criticality we therefore expect 
p to be in the range 



1 < fj, < n. 



(4) 
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Solving the stationary mean field equation = ap'"—p" 
we obtain the stationary density p = a, hence the critical 
point is a c = and the density exponent is (3 MF — 1. At 
the mean field critical point the full Langevin equation 
should be invariant under the scaling transformation 



Ax , t^A z t. 



7>- 



(5) 



where A is a dilatation factor while z = vh/v± and 
X = @/v\_ are quotients of the three standard critical 
exponents. Comparing all terms except for the noise, 
scaling invariance implies that z = 2 and \ = 2/n, i.e., 



FIG. 2: The density of particles, pit), times t s (S = 0.32) 
as a function of time for p = 0.6855, 0.6853, 0.6851, 0.6849, 
and 0.6847 from top to bottom, averaged over 1500 runs on a 
system with 4096 sites. The best straight line is obtained for 
S = 0.32 and p = 0.6851. 



Moreover, we can check the relevance of the noise term, 
which is responsible for fluctuation effects. By simple 
power counting we find that the noise is relevant below 
the upper critical dimension 



2/i 



/3 MF = 1, 



,MF 



= n/2, 



,MF 



(0) 
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while it is irrelevant above d c where the mean field ex- 
ponents (|(|) are expected to become exact. For DP 
(n = p = 1) we obtain the well-known result d c = 4, 
while for the PCPD the upper critical dimension has to 
be in the range 2 < d c < 3. This result is in agree- 
ment with recent numerical findings by Odor et al. |25( ] 
suggesting that d c = 2. 

As the main observation, which triggered the present 
work, we note that the upper critical dimension for third- 
order processes {n = 3, m < n, 1 < p < 3) is larger 
than 4/3. Consequently, in 1+1 dimensions such a triplet 
process (TP) should still be characterized by non-trivial 
fluctuation effects. Moreover, the density in the inactive 
phase is known to decay as pit) ~ (lni/i) 1 / 2 . This type 
of decay in the absorbing phase differs from all other 
previously known universality classes of phase transitions 
into absorbing states, suggesting that also the transition 
itself should belong to yet another universality class. 

b. Numerical simulations. In order to study the 
triplet process numerically, we introduce a fermionic vari- 
ant of the reaction-diffusion process (0) with n = 3 and 
to = 2. It evolves by random-sequential updates and is 
defined by the following dynamic rules: 



A -> A<D with rate (1 - p)/2 

A0^0A (l-Jj)/2 

AAA — > AA$ (l-p)/2 

AAA-+QAA (l-P)/2 

AAA% -> AAAA p/2 

%AAA -» AAAA p/2 



(8) 
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A typicial space-time plot of the process at criticality is 
shown in Fig. [l] As can be seen the process generates 
spatiotemporal structures, possibly indicating the pres- 
ence of fluctuation effects. Note that we tuned the rates 
for diffusion and particle removal in the same way as in 
Ref. @. 

Performing standard Monte Carlo simulations (see, 
e.g., 0]) we find clear evidence for a continuous phase 
transition between an active phase, where the denstity 
of particles is asymptotically constant, and an inactive 
phase, where the particle density decays algebraically 
with logarithmic corrections. Assuming that the criti- 
cal behavior at the transition obeys simple power law 
scaling we find the critical threshold p c — 0.6851(4)(see 
Fig. ||). As in the PCPD, there are strong corrections so 
that the scaling regime, if existent at all, is not reached 
before 10 4 time steps. Averaging over many independent 
runs in the time interval 10 4 < t < 10 6 we estimate the 
critical exponents by 

!/||= 2.5(2), z = 1.75(10), 6 = f3/v\\ = 0.32(1). (9) 

Similar exponents were obtained in other variants of the 
triplet process with m < n = 3 (not reported here). In 
all cases the estimates differ from the mean field (MF) 
exponents uf F = 3, z MF = 2, and S MF = 1/3, leading 
us to the conclusion that critical behavior of the 1+1- 
dimensional TP is indeed characterized by non-trivial 
fluctuation effects. As expected, these deviations are 
quite small (less than 20%) since the simulations are car- 



ried out close to the upper critical dimension. 

We note that our results are not accurate enough to 
doubtless confirm the validity of power-law scaling over a 
large range. As our simulations seem to reach the scaling 
regime only after 10 4 time steps, the accuary of our esti- 
mate for p c is limited as well. In addition, the asymptotic 
power-law behavior may be shadowed by logarithmic cor- 
rections which are also present in the inactive phase. Fi- 
nally, as in the case of the PCPD, the assupmtion of 
simple power-law scaling and the concept of universality 
may be questioned as a whole. Nevertheless we believe 
that the MF arguments and the numerical evidence are 
strong enough to conclude that this model exhibits a dif- 
ferent type of critical behavior, where fluctuation effects 
are likely to play an important role. 

Regarding the limited accuracy of numerical simula- 
tions a major drawback could be the definition of the 
model as a fermionic reaction-diffusion process with four- 
site interactions. A bosonic variant with two-site inter- 
actions is currently under investigation. Moreover, it is 
important to determine the exponent /i in the noise cor- 
relator. Preliminary simulations suggest a value close to 
/j«2. Finally, the influence of the diffusion rate has to 
be studied systematically. 
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